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problem often encountered in parking 
giirugc slabs in northern climates where 
deieing salts are heavily used. A ran- 
dom field model of the delamination 
process is investigated and compared 
against some field data. Delaniinatcd 
regions of the slab arc modeled as cx- 
cursioni of a random field above a pre- 
scribed threshold and the growth of 
these regions with time is obtained liy 
allowing the threshold to fall a.l a func- 



tion of time. .Simulation based excur- 
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mean and variability of various a.spccts 
of the delamination process using this 
model. 
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1. Introduction 

The delamination antJ spalling of concrete sur- 
faces in parking structures in northern climates is 
an ongoing and expensive problem. The Trans- 
portation Research Board [1] estimates that be- 
tween 50 and 150 parking structures in the 
Northeast and Midwest United Slates will need to 
be rehabilitated each year for the next 10 years at 
an average cost of $1 million per structure. With- 
out proper design and/or maintenance, dcicing 
salts brought in by vehicles from the roadways are 
deposited on the concrete surface along with water. 
Chloride ions gradually penetrate the concrete and 
electro-chemical processes lead to corrosion of the 
reinforcing steel. This results in both a degradation 
of structural integrity and, since the corrosion 
products occupy considerably more volume than 
the original steel, delamination or spalling of the 
concrete surface and loss of utility. 

This paper is aimed at developing a tool to aid 
in a rational probabilistic approach to the rehabil- 
itation of parking structures. Such an approach 
allows the optimal allocation of limited resources 



to this ongoing and rather expensive maintenance 
issue. A simple stochastic model involving only a 
few parameters Is used herein to represent what is 
known to be a complex phenomenon. The follow- 
ing factors are suggested by Public Works Canada 
[2J to have the largest effect on the onset of con- 
crete delamination; 

■ chloride ion input: quantities of deicing salts 
used, 

■ concrete penneabiliiy: influenced in turn by 
water/cement ratio, intensity and frequency of 
cracking, surface coatings/sealers and construc- 
tion practices, 

• ambient ternperature, humidity, precipitation 

• concrete cover depth 

• pH of local aggregates 

• conductivity: wetness of concrete. 

It is immediately recognized that virtually all of 
these factors are highly variable from structure to 
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structure and within a single structure from point to 
point. As well, on a practical basis for an existing 
structure, some of the factors are unknowable ex- 
cept through extensive destructive testing. Clearly a 
representative model should not require experi- 
mental validation on a per structure basis, where 
the experimentation may be more expensive than 
the final repairs. This, in fact, is the primary motiva- 
tion for the use of stochastic models. 

In a different approach to the same problem, 
Attwood et al. [3] develop a limit state function for 
parking structures based on a critical fraction of 
floor area delamination. Da, 



nt=Da-aS(t-t\), t>ti. 



(1) 



where m<0 denotes the "failure" state, 5 is the 
annual delamination rate, a is a correction factor 
accounting for ambient temperature and crack 
widths, / is time in years, and ii is the time to initia- 
tion of delamination (also in years). Note that the 
second term is taken to be zero for all r</|. 
Attwood et al. employ a First-Order Reliability 
Method (FORM) to estimate reliabilities associ- 
ated with the delamination of a parking garage 
structure. In their approach, all of the factors ap- 
pearing in Eq. (1) are expressed in terms of random 
variables having assumed distributions and the joint 
cumulative probabilities are evaluated using 
FORM. 

In this paper a 2-dimensional random field model 
representing the spatial delamination process over 
time is investigated. The random field may be 
loosely interpreted as the out-of-plane stress field 
at the reinforcement level (which changes randomly 
from point to point over the area of the slab). De- 
laminated regions are represented by excursions of 
the random field above some threshold which can 
be thought of as a critical concrete resistance to 
horizontal splitting. In fact the stress field itself will 
not be considered in that it is unmeasurable in prac- 
tice. Rather, the field excursions will be used to rep- 
resent the delaminated regions directly, since these 
regions are measurable to some extent. A simple 
random field model is adopted whose primary moti- 
vation is to attempt to shed light on the following 
questions; 

1) what is the mean and variability of the total 
area of delamination as a function of time? 

2) what is the average size of individual delami- 
nated regions as a function of time? 

3) how many delaminated regions can be 
expected in a slab? 



The.se are essentially questions regarding the statis- 
tics of excursion regions and so some simulation 
results regarding excursion statistics will be pre- 
sented in the next section. 

Figure 1 illustrates an example of excursions of a 
random field above some predefined threshold. In 
the context of this paper, the dark regions can be 
viewed as areas where concrete delamination has 
occurred at some fixed time. As time progresses, 
delaminated regions are expected to grow in size, 
corresponding to a falling threshold level. A falling 
threshold is equivalent to a rising mean, in the case 
of a homogeneous field. Since the excursion statis- 
tics in the next section are developed as functions of 
the threshold, the falling threshold interpretation is 
used here rather than a rising mean. In either case, 
a nonhomogeneous field could be employed, if the 
data so indicated, by considerably extending the 
simulation based study (in the absence of analytical 
results). In this preliminary investigation, only a ho- 
mogeneous field is considered. 




Fifi. I. Excursions of a Iwo-dimcnsional random field above 
threshold bu. 

To represent the delamination process, excur- 
sions of an isotropic Gaussian random field having 
mean zero and unit variance will be used. Although 
other values of mean and variance are possible, the 
excursion statistics are dependent purely on the 
distance between the mean and the threshold. This 
distance can be expressed in units of rr, the stan- 
dard deviation, so that there is no advantage in 
choosing anything other than a mean zero, unit 
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variance field. The quantity ha will be referred to 
henceforth as the physical threshold and b alone as 
the threshold. The choice of an isotropic Gaussian 
process has been made largely for simplicity, there 
being little evidence available to clearly justify other 
types of random functions. 

2. Excursion Statistics in Two Dimen- 
sions 

With respect to excursion statistics, such as the 
mean number and area of isolated excursions, ana- 
lytical results developed to date are asymptotic in 
nature, accurate only at very high thresholds where 
the excursion process approaches a Poisson point 
process. Often in engineering problems the interest 
is in thresholds which are quite a bit lower, such as 
the delamination process considered herein. In or- 
der for the proposed random field model to be use- 
ful in this context, excursion statistics should be 
available. This section summarizes a study in which 
excursion statistics are obtained through Monte 
Carlo simulation. Specifically, an ensemble of 20OO 
realizations of an isotropic zero mean, unit variance 
Gaussian random field, Z(x), with Markovian co- 
variance function, 



B(T) = o^p(T) = a^exp(-||T|). 



(2) 



are produced using the Local Average Subdivision 
(LAS) method [4, 5], where t is the lag vector and 
9 is the scale of fluctuation. The scale of fluctuation 
is loosely interpreted as the distance over which 
correlation is significant. Since many of the statis- 
tics of interest depend strongly on the scale of fluc- 
tuation, 2000 realizations were generated at each of 
5 different scales of fluctuation. 

Individual realizations are decomposed into 
excursion regions and "holes," using a space-filling 
algorithm, over a range of thresholds b =[-4, 4]. 
The mean and variance statistics of the excursions 
are estimated over the ensemble. It should again be 
eiriphasized that b is measured in units of standard 
deviation, for example b=2 implies a threshold at 
two standard deviations from the mean. For a unit 
variance field there is no distinction between the 
value ofb and the physical threshold level. However 
the local averaging performed by the LAS method 
results in a slight decrease of the variance of the dis- 
cretized field, as dictated by local averaging theory. 
In that each field is represented as a discrete lattice 
of 128 X 128 "cells," the variance of each cell varies 
from 0.971 at the smallest scale of fluctuation con 
sidered to 0,999 at the largest. The distinction 



between cell variance and point variance will be ig- 
nored in this paper, although the results presented 
in the various plots to follow are accurate in this re- 
spect. For a more rigorous treatment of this issue, 
see Ref. [6]. 

Within a given domain J^= [0, Li] X [0, Lt] of area 
Ar=L] L2, the total excursion area per unit area, 
Db, where the process Z(j) exceeds some 
threshold, can be defined by 



^""il HZ(.x)-b(T)dx, 



(3) 



where /xr is the physical threshold of interest, tr^ be- 
ing the variance of the process, and /(■) is the indi- 
cator function defined on V (taken to be zero 
outside the domain V) 



'^'^ I ifr<0 



(4) 



For a homogeneous process, the expected value of 
Db is simply 

E[D,] = P[Z^h(r], (5) 

which, for a zero-mean Gaussian process yields 

£[D,]=l-0(fc). (6) 

where 'P is the standard normal distribution func- 
tion. The estimate of £[D/i], denoted ntok^ derived 
using the simulation results is shown in Fig. 2 and is 
in complete agreement with Eq. (6). Figure 2 (b) 
shows the estimated standard deviation of Dt de- 
noted sdi,. Note that while EiDb] is independent of 
0, its variance is not. In keeping with the practice of 
normalizing all results, the scale of fluctuation has 
been normalized with respect to L =\/A^ in these 
plots. Note also that the horizontal threshold axis 
values decrease to the right— this is because the 
threshold axis is associated with time in the next 
section and time increases to the right as usual. 

Figure 3 shows the estimated mean and standard 
deviation of the number of isolated excursion re- 
gions, Nb, denoted m\^ and sn^ respectively, as a 
function of scale and threshold, Although not read- 
ily apparent from the plot, the limiting value of m^i, 
at the right edge of the plot (/>-*- -» ) is 1. In other 
words, a single excursion exists over the entire do- 
main for very low thresholds. At the other extreme 
(b-^00), the limiting value is zero as expected. 
Although the roughness of the estimate .Tatj, as seen 
in Fig. 3(b), is as yet unexplained, it appears un- 
likely that it arises from statistical uncertainty in the 
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Fig. 2. Eslimalcd stalistic^ uf total excursion aiea per unit area, Db,: a) mciin, b) standard dcvialion. 
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Fig. 3. Estimated stalislin of number of isolaltd excursions, N^: a) mean, b) standard dcvialion. 



estimation procedure. Based on a sample size of 
2000 realizations, a 90 percent confidence bound on 
the standard deviation of TVs is about ±0.025s,vj, 
which is smaller than the size of most of the 
"bumps" seen in Fig, 3(b). 

Within a given realization, the average area of an 
isolated excursion per unit domain area, Dc, can be 
obtained using the number of excursions. 

Since Db is the sum of the Nt isolated excursion 
areas, the expected value of De is just 

Figure 4(a) illustrates this result using the esti- 
mated mean value of M shown in Fig. 3(a). The 
estimated standard deviation of Dc shown in Fig. 



4(b) is derived using the assumption that the sizes 
of isolated e.\cursions are independent. While this is 
true from realization to realization, it is clearly not 
true within a single realization. Thus Fig. 4(b) can 
only be considered to be a rough indication of the 
true variability of the area of isolated excursions, 

3. Calibrating the Model 

There are essentially only two parameters in the 
random field model considered in this study. These 
are the scale of fluctuation, 6, and the threshold 
level, h. Ahhough there is little published experi- 
mental evidence to allow a clear statement of what 
these parameters should be for a given structure, 
some preliminary estimates are possible. First the 
relationship between the threshold level and the 
age of the structure can be obtained by fitting 
delamination versus age data collected by Trow [7] 
and presented in Ref. [3] to the curve of Fig. 2(a). 
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A good fit was obtained using least squares re- 
gression (r^ = 0.93) on the linear relationship 
b= 3.17— 0.225r. To obtain this relationship, the 
total fraction of delaminated area observed by 
Trow, Db, was plotted against h = '(1 ~Db) and 
a best fit line obtained. Figure 5(a) shows the re- 
gression results and Fig. 5(b) shows where the ob- 
served delamination results would appear on Fig. 2. 
To expand the scale, only positive thresholds are 
shown in Fig, 5(b). The choice of function relating 
time and threshold is largely arbitrary as long as 
b(t) is a decreasing function over all times of inter- 
est (assuming delaminated areas cannot "heal"). 
While a quadratic gives a slightly better fit to the 
raw data, it violates this principle and so cannot be 
used to extrapolate. 



Some additional unpublished data was made 
available to the author by Public Works Canada 
based on a survey of a single parking garage struc- 
ture in Ontario, Canada. Figure 6 illustrates this 
data. The ±1 standard deviation curves are ob- 
tained using Fig. 2(b) with 6/L=0A. Clearly, for 
such a choice in 0, the variability in Db predicted by 
the random field excursion model underestimates 
the variability in the observations. The additional 
variance arises because the time-threshold relation- 
ship b{t) is itself a function of random coefficients, 
as implied by the regression analysis. Alternatively 
and equivalently, the additional variability could be 
ascribed to the fact that the delamination field is 
not homogeneous on such a scale. When considered 
at the scale of a typical bay, some bays show much 
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Fig. 4. E&timatcd statistics of isolated excursion areas per unit area, Dc,: a) mean, b) standard deviation. 
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Fig. J. Regression results using Trow's delamination rate data: a) regression, b) observed versus predicted D*. 
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higher delamination rates than others, perhaps cor- 
responding to smaller mean cover depths and/or 
higher Cl~ inputs. Note also that larger values of 61 
L lead to higher standard deviations in Dt, so that 
the variability seen in Fig. 6 could also be explained 
by larger values of 0. However the choice of 
should not be based on the variability in Db, rather 
this observation indicates that perhaps the Gauss- 
Markov covariance model is not appropriate. 

In this study, ^(r) will be treated as a determinis- 
tic function which corresponds to the choice of a 
homogeneous random field. It is beyond the scope 
of this initial investigation to consider estimating 
the parameters of a nonhomogeneous (or self- 
simitar) field to represent the delamination process 



generally, although this appears indicated. Since 
plan views of the delaminations are available from 
Public Works for individual bays of approximately 
5m X 5 m in extent (L =5), attention will be re- 
stricted to data on such a scale over which the field 
can be considered homogeneous. Choosing one 
such area from the Public Works data, the mea- 
sured delamination fraction as a function of time is 
shown in Fig. 7. In this case the line of best fit was 
found to be 



fc=4.13-0.172f 



(V) 



with r^ = 0.86. Notice that for such a case, the vari- 
ability in the observations is substantially reduced 
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Fig. 6. Delamination rate data provided by Public Worki Canada. The fitted line is ij =3,61 -^0.142/ with /■^ = 0.6. 
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Fig. 7. Deliminatinn rate provided by Public Works Canada over a single slab region. 
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and that 2/3 of the observations lie on or within the 
±1 standard deviation curves for 0/L=0.] (see 
also Fig. 5). This result appears encouraging al- 
though it is recognized that it could be due in part 
to the reduced number of samples, even though 
similar results were found for most other bays. 

Before considering the estimation of fl, it is worth 
pointing out a further difference between the ran- 
dom excursion model and the commonly accepted 
delamination model. Salt induced reinforcement 
corrosion is generally believed to involve two stages 
(see Ref. [8]): 1) an initiation phase during which 
the alkalinity of the concrete surrounding the rein- 
forcement {which renders the steel passive) is re- 
duced by the migrating 0~ ions, and 2) an active 
phase during which corrosion takes place. Attwood 
et al. [3] estimates the initial phase to last 4.7 years. 
During this phase no corrosion is assumed to take 
place. Note that such a model can only be applied 
to points in the slab where the concrete is in contact 
with the reinforcement. 

In contrast with this two stage model, the random 
field excursion model admits some probability of 
delamination even at time f =0. Using Eq. (7), one 
obtains h =4,\3 at time r =0, so that the expected 
total delamination area per unit area is 1.7 x 10"^ or 
about 17 mm^/m- of slab. At these levels the precise 
definition of delamination comes into question. If it 
is strictly interpreted as a loss of bond between the 
concrete and reinforcement then this result is not 
unreasonable given the presence of cracks, voids 
and the initial state of the reinforcement. At time 
/ =5 years, the expected total delamination area is 
still only about 6 cmVm^ of slab, a level which is 
probably still largely undetectable at the surface of 
the slab and presumably would correspond to corro- 
sion in the immediate neighborhood of surface 
cracks. Although it is believed that cracks do not 
contribute significantly to the areal delamination 
process [8], it is not unreasonable to expect that 
they can be initiators of the corrosion process at dis- 
crete points in the slab. If this is the case, then the 
many year delay before the onset of observable lev- 
els of delamination implies that the corrosion 
growth should be quite slow at first. The results 
predicted by Fig. 2(a) are in basic agreement with 
this in that £[/>*} grows very slowly for b decreasing 
to about 2 [corresponding to f < 12 years using 
Eq. (7)]. 

Turning now to the estimation of the scale of 
fluctuation, 6, it becomes apparent that this task is 
complicated by the type and quality of data avail- 
able. Ideally, one would take measurements of the 
corrosion induced stress field over a number of 



structures, estimate a spatial covariance structure 
and from this obtain 0. Even if this approach were 
possible, the nonhomogenieties mentioned above 
would make it difficult. However, in general the 
stress field is unmeasurable and what little data is 
available generally consists of surveys giving the 
spatial extent of delaminated regions. Figure 8(a) il- 
lustrates such a survey while Fig. 8(b) is a realiza- 
tion of the random field excursions using 6 = 0.5 m. 
Once the time-threshold relationship has been 
established, and for the purposes of this argument 
Eq. (7) will be used, a possible technique of estimat- 
ing 8 would be to count the average number of 
excursions and enter Fig. 3(a) at the appropriate 
threshold to estimate 6/L. Purely on the basis of 
Fig. 8(a) this yields estimates of = 2 to 4 m (0/ 
L-OA to 0.8). However realizations at this scale 
yields excursions which are generally far too well 
connected as shown in Fig. 9. Realizations at such 
large scales appear like large land masses with many 
small islands close to shore. Figure 8(a) has 
"islands" that are more uniform in size and distri- 
bution implying a smaller scale of fluctuation. In 
Fig. 8(b), produced using 8 = 0.5, the larger islands 
are of similar size to the dclaminations seen in Fig. 
8(a). This along with arguments to follow supports 
the choice of a smaller scale of fluctuation. 

While one cannot expect the plots in Fig. 8 to be 
identical since they are both independent realiza- 
tions of a random process, any more than one could 
expect the pattern of dclaminations in another 
building to be identical, a number of points can be 
made about the two plots; 

1) there is no apparent spatial orientation of the 
delamination regions in the observations of Fig. 
8(a), indicating that the assumption of isolropy is 
acceptable, at least for this case. 

2) about half of the excursions in the random field 
model (Fig. 8b) are of very small extent. This frac- 
tion increases at larger scales. On the other hand, 
in Fig. 8(a) there are only very few "small' areas 
appearing in the later survey. It seems reasonable to 
suspect that additional small delamination regions 
are in fact occurring in the real slab but that the 
chain-drag surveying technique is unable to resolve 
them. Operator bias will almost certainty also be 
present due to the prior knowledge of existing de- 
lamination areas. 

3) the random field model is much "rougher" 
than the observed delamination plots. Again this is 
likely a problem with the ability of the chain-drag 
survey to resolve detail. 
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Fig. 8. a) Observed dclafnin;ilitin regions at ages ( = 13 years (dark grey) and ; =22 years (lighl grey) on a 
5mx5ni portion of a parting garage slab, b) Excursions of a random field above i>-0J4(/ = 13) and 
ii = 19(/-22) using 6=0.5 and Eq. (7). 




Fig. 9. A realization of excursions at f » 13 years (dark grey) 
and at/ =22years (light grey) using = 4 mim a 5 m x 5m ficlil. 



The last two points illustrates the difficulty in esti- 
mating 6 on the basis of delamination surveys. The 
chain-drag method depends on setting up reverber- 
ations in the delaminated concrete. For delamina- 
tion details below a certain size, the frequency .shifts 
are undetectable to even the most sensitive human 
ears, rendering them unnoticeable. Thus many of 
the delamination details on which an estimation of 
6 depend are unavailable using current surveying 



techniques. To some extent the deficiencies in the 
chain-drag procedure could be accommodated sim- 
ply by introducing more local averaging in the ran- 
dom field model — effectively smoothing the field. 
While such a "correction" would not Hkely result in 
an improved delamination model, it may allow 
improved estimates of 0. For the purposes of this 
investigation, the realization based estimate of 
^ = 0.5 m is used. 

4. Discussion 

Conceptually the model proposed herein is quite 
attractive in that the delamination process is indeed 
a threshold excursion process in two dimensions. 
Once the details of the model have been established 
(type of distribution, time-threshold relationship 
and scale of fluctuation) and some of the properties 
of threshold excursions in two dimensions have 
been determined analytically or via simulation, the 
model can be used in a reliability context. For ex- 
ample, using Eq. (7), at time / = 12 years {b =2.07) 
Fig. 2 along with Eq, (6) indicates that widj = 0.0194 
and5/3;^ =0.0080, using fl = 0.5 m. If D^ is assumed to 
have a Beta distribution, which is properly bounded 
between and 1, then the probability that Dt, is less 
than 3 % at 4 years is given by 



/'[Oi,<0.03] = 0.90 



(8) 



In turn, if a target reliability of 0.9 were chosen 
against delamination in excess of 3 %, then 
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inspection of the garage would be recommended at 
r = 12 years for this structure. 

Because the model includes the spatial aspects of 
delamination, it can be used to evaluate testing and 
mapping procedures. For example, at time t = 12 
years one could expect about 42 isolated excursions 
on a 5 m x5 m slab (0 =0.5 m) with a standard de- 
viation of about 11. The area of each isolated excur- 
sion would average about 0.011 m^ with a standard 
deviation of about 0.022 m^. Sun'eys which yield re- 
sults considerably different than these may require 
verification using other techniques. In addition, the 
spatial description of delamination could be used in 
a structural reliability study. For this the simulation 
approach could be employed to yield a measure of 
the degree of clustering of the delamination regions 
(see Ref. [6]). 

One recognizes that the excursion model is at- 
tempting to predict the cracked state of a concrete 
slab. In that internal cracks are exceedingly difficult 
to map, even in controlled laboratory conditions, 
the model is to some extent intuitive and will likely 
remain so until improved surveying techniques are 
developed. Nevertheless the model demonstrates 
some promising features and can be used as a pow- 
erful reliability tool when its parameters are clearly 
defined in terms of additional data and simulation 
studies. In particular, the data shown in Fig. 6 indi- 
cates that perhaps alternative correlation functions 
[see Eq. (2)] should be studied — multiple-scale or 
self-similar type random fields are suggested, show- 
ing the small scale behaviour over small regions 
while reflecting also the large scale, slower varia- 
tions over larger domains. Also the fact that the 
model allows delamination to occur at time ;= (or 
before) implies that some thought should be given 
to the assumption of a Gaussian random field and/ 
or the time-threshold relationship. In the interim, 
however, the choice of a Gaussian field and linear 
time-threshold relationship leads to results which 
appear reasonable for any time r >0. 
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